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Abstract 

Dynamics of four-dimensional massless fields of all spins is formulated in the 
Siegel space of complex 4x4 symmetric matrices. It is shown that the unfolded 
equations of free massless fields, that have a form of multidimensional Schrodinger 
equations, naturally distinguish between positive- and negative-frequency solutions 
of relativistic field equations, i.e., particles and antiparticles. Multidimensional 
Riemann theta functions are shown to solve massless field equations in the Siegel 
space. We establish the correspondence between conserved higher-spin currents in 
four-dimensional Minkowski space and those in the ten-dimensional matrix space. It 
is shown that global symmetry parameters of the current in the matrix space should 
be singular to reproduce a nonzero current in Minkowski space. The I?— function 
integral evolution formulae for 4d massless fields in the Fock-Siegel space are ob- 
tained. The generalization of the proposed scheme to higher dimensions and systems 
of higher ranks is considered. 
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1 Introduction 



The idea that a set of massless fields of all spins in the four- dimensional space should 
admit a natural description in the ten-dimensional space with real symmetric matrix 
coordinates X^^ = X^^ {A, B, . . . = 1, . . . ,4) was originally proposed by Fronsdal in [1]. 
Later, the same conclusion was independently reached in |2]. The dynamical equations in 
J^M, that for M = 4 are equivalent to the field equations for massless fields of all spins 
in the four-dimensional Minkowski space M^, are very simple [3]. All integer spin fields of 
are described in A^4 by a single scalar field b{X) that satisfies the Klein-Gordon-like 
equation 

All half-integer spin fields are described by a single fermion field /b(^) that satisfies the 
Dirac-like equation 



^Jc{X)-—-^fB{X) = Q. (1.2) 



The equations fll.lj) and (11.21) were derived in |Sj from the system of equations 

d 92 



ax^i.+/^^yI^J^(^W = 0' (1.3) 

where were treated as auxiliary commuting variables (the parameter /i 7^ is in- 
troduced for the future convenience). Although the equations (ll.ip -f fTTSl) were originally 
considered for M = 4, they make sense for any M. 

The equations (11. 3p express the first derivatives with respect to space-time variables 
X^^ in terms of the fields themselves. As such, they belong to the class of unfolded partial 
differential equations (PDE) that, more generally, express the exterior differential of a set 
of differential forms in terms of exterior products of the differential forms themselves. Such 
a first-order form of dynamical field equations can always be achieved by introducing a 
(may be infinite) set of auxiliary fields which parameterize all combinations of derivatives 
of the dynamical fields that remain non-zero on the field equations. For example, in the 
system (11.31) . the dynamical fields are 

h{X) = C(0|X) (1.4) 

and 

/AW = ^^m^)L- (1-5) 

As a consequence of (11.31) . they satisfy, respectively, the equations (II. ip and (II. 2p . All the 
fields 

gA.....„(X) = ^^,^^^^^^,^ C(y|X)|^^^, n>l (1.6) 

are auxiliary, being expressed via higher X-derivatives of the dynamical fields by virtue 
of the equations (11.30 . In [3] it was shown that the equations (II. ip and (II. 2p along with 
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constraints that express the auxihary fields via X-derivatives of the dynamical fields 
exhaust the content of the unfolded system {\1.3\i . That the equations (11.31) formulated in 
the ten-dimensional space-time, still describe massless fields in four dimensions was also 
shown in [3] using the unfolded dynamics approach (see also Section [2]). 

Theories in A4m have been studied in a number of papers from different perspectives 
[1[51[61[71[81[91[T01[II1[I21[I31[T1[I51[I61[I3 In this paper, we will further study 

the higher-spin (HS) theory in the matrix space. 

The main practical goal is to show how HS conserved currents in A^4 found in [5] 
reproduce usual HS conserved currents in Minkowski space found in [20]. The analysis is 
not completely trivial since the conserved charges in contain an additional integration 
over one spinning variable. The apparent difficulty is that the compact spin space is 
contractible to zero, hence implying that the charge must vanish for regular solutions of the 
field equations in A^4. A standard way out would be to integrate over a noncontractible 
cycle in AiM- (In fact, in its Sp{2M) invariant compactification which is Lagrangian 
Grassmannian PQ S].) We have not been able to proceed along these lines, however. 
Instead we will show in this paper that the Minkowski charge is correctly reproduced by 
virtue of introducing a singularity that effectively makes the integration cycle over the 
spinning variable noncontractible. The obtained results may have several applications. 

First of all, conserved currents determine the lowest order Noether interactions with 
the HS gauge fields associated with the HS symmetries. It is straightforward to introduce 
cubic interactions of HS gauge potentials with conserved currents via replacing the global 
symmetry parameters rj in the charge by the corresponding HS one-form gauge connec- 
tions. The results of this paper show that, to reproduce correctly the HS interactions in 
the four-dimensional setup, HS potentials in should develop a singularity in the spin- 
ning directions. In other words, the results of this paper indicate that there are nontrivial 
fiuxes in the spinning directions in AI4, that support charges in Minkowski space. 

Another application is that the obtained formula for conserved charges allows us to 
write the integral evolution representation for solutions of massless field equations in A44 
with the help of D-functions introduced in [4j. Generically, P-functions provide the 
integral representation for solution of field equations of the form 



where E is a surface where the initial data are given. Formulae of this type should respect 
a number of properties. Firstly, V{X,X' ) should form a solution of the field equations 



the formula (11.71) should be independent of local variations of S which property is satisfied 
if T>{X, X' ) solves the field equations with respect to X' and (I1.7P is defined as a conserved 
charge with respect to X' . The proper definition of the integration measure dX' in (11.71) is 
achieved in this paper. Note that so defined P-function satisfies the composition property 




(1.7) 



under consideration with respect to X. Secondly, V^X, X') 



X,X'&S 



(5s(X-X'). Thirdly, 
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In the analysis of HS currents, we find it most convenient to depart from the real 
space M.M to its complexification S^m known as upper Siegel half-space [21] (see also 
[22j). It turns out that in this framework positive- and negative-frequency solutions 
identify with holomorphic and antiholomorphic solutions in the upper Siegel half-space. 
J^M is a boundary (absolute) of the Siegel half-space. A surprising conclusion will be that 
unfolded field equations themselves distinguish between positive- and negative-frequency 
solutions of the field equations, i.e., between particles and antiparticles, the property 
usually delegated to the quantization procedure. It is worth to mention that the unfolded 
equations have a form of multidimensional Schrodinger equations. These conclusions may 
eventually be of key importance for a deeper understanding of the interplay between 
unfolding and quantization. 

Once the HS field equations are reformulated in the Siegel space, it is straightforward to 
observe that they are solved by Riemann theta functions. This fact is truly remarkable in 
view of the role that theta functions play in modern geometry, theory of integrable systems 
and String Theory, and is hoped to shed more light on fundamental structures underlying 
HS theory. Let us note that in some sense theta functions provide most symmetric non- 
zero solutions of massless field equations. Namely, they are invariant up to a phase under 
the transformations from the Igusa group ri^2 C Sp{2M,Z) (see [22] for more detail). 
This class of solutions may indeed play a distinguished role in the HS theory because the 
observables constructed from such solutions, like, e.g., conserved currents, turn out to 
be invariant under Fi 2. Note also that some of the properties of theta functions admit 
natural interpretation in terms of the unfolded massless field equations. 

All seemingly different aspects of HS theory considered in this paper take their origin 
in the symmetry properties of the massless field equations (11.31) . i.e., HS symmetries. One 
of the advantages of the unfolded formulation is just that it makes symmetries of PDE 
manifest. Therefore we start in Section [2] by recalling some relevant facts of the unfolded 
formulation approach, with the emphasize on symmetries in Subsection 12. 1[ In Subsec- 
tions 12.21 and 12.31 we recall, respectively, the relationship between the ten-dimensional 
space Ai^ and Minkowski space M*^ and some known results on HS conserved currents in 
A^4 and M*^. In Section [31 we extend the unfolded description of massless higher spins 
to the complex Siegel space and explain how unfolded field equations distinguish between 
positive- and negative-frequency solutions. In Section HJ we develop further the construc- 
tion of conserved currents by extending Ai m to Mm x x I^*^ where we introduce a 
2M-form, that is closed by virtue of certain unfolded equations, and show how it repro- 
duces the previously known closed M-forms associated with conserved currents. Then, in 
Section, [S] we show how the known 4d currents result from those in 7M4 after introducing 
a singular fiux in the spinning directions. Using P-functions found in pE] and the obtained 
construction of conserved currents, we derive in Section [6] the integral evolution formulae 
for solutions of the massless field equations in Ai^. Multilinear conserved currents are 
considered in Section [71 In Section [SI we show that Riemann theta functions form a nat- 
ural class of periodic solutions of massless field equations. Conclusions and perspectives 
are discussed in Section [91 In Appendix, we describe a commutative associative product 
law o which endows the space of solutions of unfolded equations of any rank in Mm with 
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the commutative ring structure. 



2 Preliminaries 

A natural approach to the study of dynamical equations of motion in the HS gauge 
theory, referred to as unfolded formulation, consists of reformulation of PDE in the form 
of certain covariant constancy conditions [23]. Using this approach, consistent gauge 
invariant nonlinear HS equations of motion were found in [2H [251 [26] for HS theories 
in three, four and any dimension, respectively, (see [271 [28[ [29] for reviews and more 
references). The unfolded formulation is particularly useful for revealing symmetries and 
dynamical content of PDE as discussed e.g. in [10]. Here we briefly recall some properties 
of this approach. 

2.1 Unfolded formulation and symmetries 

Consider a system of linear PDE of the form 

{d + uj)C{X)=0, d = dX''^, (2.1) 

where C{X) is a section of the trivial vector bundle B = M.'^ xV over the space-time base 
M'^ with the local coordinates X'' 

V — > B 
i 

with a linear space V as the fiber. In the cases of interest V identifies with an appropriate 
space of power series f(Y) = X^^o fAi...A„{^)y^^ ■ ■ ■ Y^" in some auxiliary variables F^, 
i.e., C{X) with values in V is realized as a function C{Y\X) of the two types of variables. 
The one-form OJ^X) = dX^uj^^X), that satisfies the flatness condition 

rfcu + ^[cu,cj] = 0, (2.2) 

is some fiat connection of a Lie algebra g C End V with the Lie product [ , ] . (Here we 
discard the wedge product symbol A .) 

The equation (12. ip is invariant under the global symmetry q. Indeed, the system (12. ip 
and (12. 2p is invariant under the infinitesimal gauge transformations 

5uj{X)=dt{X) + [uo{X),t{X)], 

5C{X) = -e{X)C{X) , (2.3) 

where €{X) is an arbitrary symmetry parameter that takes values in g. For a given uj{X), 
there is a leftover symmetry with the parameter e(X), that satisfies 

6uj{X) = de{X) + [uj{X), e{X)] = . (2.4) 
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This equation on e{X) is consistent as a consequence of fl2.2p . i.e., the Bianchi identity 
d'^ = does not impose any further conditions on e(X). Therefore, it reconstructs locally 
the dependence of e{X) on X in terms of its values e(Xo) at any point Xq of space-time. 
In the absence of topological obstructions, the resulting global symmetry algebra with 
the parameters e(Xo) is q. It is therefore enough to observe that some dynamical system 
can be reformulated in the form (12.11) . where a flat connection uj{X) takes values in some 
algebra g that acts in V, to reveal the global symmetry q (12. 3p of the system (12. ip . This 
approach is general since every g-invariant linear system of PDE can be reformulated in 
the form (12.11) by adding enough auxiliary variables (nonlinear systems are described in 
terms of an appropriate generalization associated with free differential algebras [301 EI] 
as explained e.g. in jH |19l [231 [27] ) . 

From this analysis it follows [3] that the equations (11.31) and, hence, (II. ip and (ll.2p are 
invariant under sp(2M, M) and its infinite-dimensional HS extension. Indeed, the 0-form 
C{Y\X) can be interpreted as a section of the fiber bundle B = R~^t^ x V where V 
is the Fock module of the associative Weyl algebra Am with the generators and Pa, 
that satisfy 

[Pa,Y'']=6^, [Y^,Y^] = 0, [Pa,Pb] = 0. 

The Weyl algebra is spanned by various polynomials a{Y, P). The Fock module V is 
spanned by the vectors 

V: f{Y)\0) 

generated from the Fock vacuum |0), that satisfies Pa\0) = 0. Clearly, V forms a mod- 
ule of Am as well as of the Lie (super) algebra (s)[)s(2M) constructed from Am via 
(anti)commutators with even and odd subspaces identified with the spaces of even and 
odd functions a{Y, P), i.e., a{—Y, —P) = {—l)'^^"-^a{Y, P) . The algebras of this type were 
identified in [32] with the HS symmetry algebras found in [33jl|. Note that we do note rule 
out the Lie algebra f)5 resulting from A^^ via commutators for all its elements because, 
although, having wrong relationship between spin and statistics, it is of interest in the 
context of the theory of theta functions considered in Section [HI Note also that bosonic 
spinorial symmetries of this class were recently considered in the four-dimensional setup 
in [m [Ml [37]. 

The finite-dimensional subalgebra sp(2M , M) C (s)f)s(2M) is spanned by the genera- 
tors 

La"" = \{Pa,Y^], Pab = PaPb, K^^ = Y^Y^. 
The equations (11.30 now take the form (12. ID with 

d = dX^''-^, u; = ^idX^''PAB- (2.5) 



^Note that later on it was shown [51] that, although the even subalgebra of the Weyl algebra is indeed 
the HS symmetry algebra of a nonlinear bosonic HS model with integer spin fields, in presence of fermions 
both the symmetry and the field spectrum have to be doubled by adding the Klein operators [35] (see 
[19j for a recent discussion in the context of the analysis in Mm)- In this paper we do not introduce the 
Klein operators, which however should be expected in a nonlinear supersymmetric theory. 
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Since [Pab,Pcd] = 0, so defined u; is a fiat sp(2M) connection, that satisfies (12.21) . 
Locally, any fiat connection admits a pure gauge representation 

^(P, Y\X) = g-\P, Y\X) dg{P, Y\X) . (2.6) 

Using this representation, we solve (12.41) in the form 

e(P, Y\X) = g-\P, Y\X)eo{P, Y)g{P, Y\X) , (2.7) 

where eo{P,Y) is an arbitrary X-independent element of Am 

eoiP, Y) = eJ2 VM...a/'-''-Y^' • • • Y^-Pb, . . . Pb^ ■ (2.8) 

n,m>0 

For the fiat connection (12. 5p . the pure gauge representation (12. 6p holds with 
g{P, Y\X) = exp (fiX^^'PAPB) , g'^P, Y\X) = exp (-/iX^^P^P^) • 

Let eo(P, Y) be of the form 

eo(P,r) = -eexp/i^Pe expjA^^ = -eexp j^/i^ expj^F^ exph^Ps, 

where Ja and h"^ are numerical parameters. Then the global symmetry parameter (12. 7p 
takes the form 

e(P,r|X) = -eexpjAh''exp{-fiX^''jAjB)expjAY^ exp((/i^ - 2/iX^^jc)PB) • 

The resulting global symmetry transformation (12.31) reads as 

6C{Y\X) = eexpjAh^ exp {jaY^ - iiX^^'jaJb) C{Y'' + - 2/iX^^jc|X) . (2.9) 

Note that this formula was derived in Section 7.2 of [3J in the Weyl ordering, i.e., for 
/X = 1 and 

e[,(P, Y) = exp(j^y^ + /i^Ps) = exp Pj, exp j^F^ exp(-ij^/i^). 

Differentiating (12. 9p with respect to and ja-, it is easy to derive the transformation 
law for any global HS symmetry with polynomial parameter eo{P,Y). In particular, at 
/i = 1 and M = 4, the sp(8) transformations are [3] 

K^''CiY\X) = ^ e-HC{Y\X) 

OJaOJb hA=jA=0 

= (y^'Y^ - 2F^X^^^ - 2F^X^^^ - 2X^^ - 4X^^X^^^^) C(r |X) . 

Note that, as explained in more detail in [3l|T9], sp(8) contains 4d conformal algebra 
su(2, 2) as a subalgebra. Therefore the equations (II. 3p are conformal invariant. Irreducible 
5u(2, 2) invariant subsystems correspond to 4d massless fields of different spins. 
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2.2 Initial data problem 

Usual d— dimensional Minkowski space-time M'^ is a subspace of the matrix space M.m 
for an appropriate M. To describe the embedding of 4d Minkowski space-time into Ai^ it 
is convenient to use complex notations with two-component indices a, (3 and a', (3' in place 
of the four- component indices A,B . . . with the convention that the complex conjugation 
interchanges unprimed indices a,P = 1,2 with the primed ones a',P' = l',2'. We use 
notation with four-component indices being equivalent to a pair of primed and unprimed 
two-component Greek indices (e.g., A = a, a') and = (Y°', ), where = Y^. 
In terms of two-component complex spinors we set 



so that X""'^' is complex conjugated to X"^, i.e., = X"'^', while X"^' is hermitian, 
j^ap' — j^fSa ^ Pqj, Minkowski coordinates, we will sometimes use notation x"^ instead of 
X°^ . Two-component indices are raised and lowered according to 

j^a _ ^O'P ^ = SapA"' , Sap = i ^12 = 1 5 

and analogously for primed indices. 

Minkowski time t and space coordinates are 

X^P' = tT'^P' + x'af , i = 1, 2, 3 , (2.10) 

where T°^' = 5°"^' while o""^ are hermitian traceless Pauli matrices. The Klein-Gordon 
and Dirac equations in Minkowski space read as 

,>(X)=0, (2.11) 



^ - ^UX) = , (X) - ^MX) = . (2.12) 



As shown in [HE], the generalized space-time Mm admits the well-defined no- 
tions of future and past. The (past) future cones of the origin X = are formed by 
(negative)positive-definite matrices X^^. There is a single time evolution parameter 

i=^X^^^AB, (2.13) 

where is some positive-definite time-arrow matrix. For a chosen time parameter t, 
the global space-like Cauchy surface Sf is parameterized as 

X^^ e St : X^^ = + tr^^ , (2.14) 

where the space coordinates x^^ are arbitrary T— traceless matrices 

x^^r^B = o, r^Br^^ = 5^. (2.15) 
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A particular solution in Ai m can be reconstructed from the values of the fields along 
with some their time derivatives on the global Cauchy surface. However, because the 
system of equations (11. ip on the field h{X) is overdetermined, some of these equations 
play a role of constraints that restrict the choice of the initial data on the global Cauchy 
surface. Independent initial data can be given on a lower-dimensional object called local 
Cauchy bundle E, which is a M— dimensional fiber bundle over a (c? — l)-dimensional base 
manifold cr G S treated as the space manifold. The Minkowski space-time is x cr C A^m 
where i? is a time axis. 

To see that initial data for the equations (11.11) . (II. 2p should be given on a M- 
dimensional surface is most convenient by using their unfolded form (11.30 . Indeed, the 
generic solution of (II. 3p can locally be given in the form 

Ciy\X) = exp (^-f^X^-^^^^ CiY\0) , 

where the "initial data" C(y|0) is an arbitrary function of M variables Y"^. 

Note that the equivalence of the unfolded equations (11.31) to (II. ip . (11.21) shown in [5] 
manifests itself in the inverse formula that reconstructs the dependence of C{Y\X) on Y 
in terms of any functions b{X) and /a^X) that satisfy (11.10 and (II. 2p . respectively. 



C{Y\X) = cos(.)6(X) + '-^Y^UX) , . = J-Y^Y^-^ . (2.16) 



2.3 Higher-spin currents 



The infinite set of conformal HS symmetries found in |3] is parameterized by various global 
symmetry parameters (12.81) . This suggests the existence of the corresponding conserved 
HS currents. Indeed, in J5j it was shown that the M-form in A^a/ 

^]M(r/,C^CO = ec,...CMrf^''^^^ A...ArfX^«^^^ (2.17) 

IBi...Bn ■■■ Ai...AM + m+n^ ''' 

where eci...CM totally antisymmetric multispinor and constants r] are the HS sym- 

metry parameters, is closed provided that the generalized stress tensor ^^{X) 

T''m...aAX) = ^ • • • g^C\Y\X)C\^\X) (2.18) 

is built of the fields C^(y|X) and C\iY\X) that satisfy ([OD- (Here k and / are color 
indices which take an arbitrary number of values.) The charge 

g(r?, C\ = j nMiv, C\ C) , (2.19) 

EM 

is independent of local variations of a M-dimensional surface E^^ , i.e., it conserves. 
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On the other hand, HS charges in Minkowski space have the form 

Q{V,C\&)= J ^]rf_l(r/,C^C'), (2.20) 

where fld-iiv^ C*^, C') is a on-mass-shell closed (ci — l)-form dual to the conserved current, 
and a'^~^ is a (rf — l)-dimensional surface in the Minkowski space-time, usually identified 
with the space surface R'^~^. The explicit expression for the on-shell closed three-form 
VL^{ri ^ , C^) in Ad Minkowski space, obtained recently in [20], is 



y=o' 



n^{r], C\ = dxaa' A dx"^' A dx^'''w^w^>7]{w, u)C''{Y\x)C\iY\x) 
where 

_ ^ _ ^ a _ aa' ^ a' _ aa' ^ 

~ ' ~ dV^' ' ^ ~ ^ Qya' ' ^ ~ ^ Q ' 

Equivalentl}|§, 

VL^{ri, C^ C') = dxr^^, A dx'^^' A dx^"^' 

ila{n)a'{m) X ...X X ...X J- ^ I3{p) ^i{m) '-/' f3' {q) ^i' {n) ) 

is closed provided that T^' restriction of the stress tensor fl2.18p to the 

Minkowski space C that has only nonzero coordinates X""^' = x"^' among X"^^, 
and the fields C''{Y\x) satisfy the id unfolded equations 

\c\Y\x) + —^C\Y\x)=0. (2.22) 



g^a/S' VI/' dY'^dY'^' 

Let us note that the equations of this type naturally appear in the study of so-called 
twistorial world-like particle models (see e.g. [381 EH]) as the Fourier transformation of 
the Dirac constraints on the respective phase space momenta. 

It was conjectured in [5] that the charge (12.191) at M = 4 should reproduce the 
Minkowski HS charge via an appropriate reduction of Ai^ to M^. This conjecture 
sounds plausible because the two charges contain the same symmetry parameters and 
the Minkowski stress tensor results from the restriction of the generalized stress tensor 
in to the Minkowski space. One of the goals of this paper is to establish the precise 
correspondence between the Ai^ and Minkowski realizations of the conserved charges. 

For the case of M = 2 the identification is a sort of trivial because = M.2 and 
the local Cauchy bundles are the same, both being two-dimensional. (Note that since 
"id conformal Minkowski charges are constructed from 3c? massless scalar and spinor they 
coincide with the 3(i conformal HS currents found in [30].) For higher dimensions the 
precise identification is less trivial. The problem is that the dimensions of local Cauchy 
bundles are different in M'^ and M.m for M > 2. The extra dimensions in compared to 



^We use the convention where a number in parentheses next to an index denotes a number of sym- 
metrized indices. For instance, a{n) stands for n symmetrized indices ai . . . an 
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(7°' ^ are responsible for spin and are associated with the corresponding compact spaces. 
It turns out |2l H [13] that E^' = M^^^ x 5^^^"'^+^ for M = 2,4,8,16 correspond to 
(i = 3, 4, 6, 10. In particular, the M = 4 spin space is S^. 

3 Quantization and Siegel space 

The coefficient fi in front of the second term in the unfolded equations (11. 3p was irrele- 
vant within an expansion in powers of Y"^ used in [3]. Its absolute value can certainly 
be normalized arbitrarily by a rescaling of Y^. However, its phase should respect real- 
ity conditions. Surprisingly, it distinguishes between positive and negative frequencies, 
i.e., particles and antiparticles. We interpret this observation as an indication that the 
unfolded dynamics encodes quantum physics. 

Indeed, general solution of the equations (11.11) and (11.21) is |1] 

biX) = 6+(X) + b-{X) = (3.1) 
^|d''^{b^Oew{^U^BX^''} + r(Oexp{-zeAeBX^^}) 

and 

fA{X) = fA+{X) + fA^iX) = (3.2) 
-^Jd"^ U{r{Oew{^U^BX^''} + /-(Oexpl-za^BX^^})- 

Note that 6±(A) = (-1)^6±(-A) and /±(A) = (-l)*'^+7±(-A). In particular, for even 
M, &^(0 and /^(O cire even and odd functions of ^, respectively. The integration in 
(13. ip and (13. 2p is hence over M.^'^ /Z2. The point = is invariant under the Z2 reflection 

— > —C,a and therefore is a singular point of the orbifold M.'^ /Z2. 

Now we observe that the unfolded equations 

{J^±'n^^)c-(Y\x) = o (3.3) 

distinguish between the positive- and negative-frequency solutions 

C^{Y\X) = -^1 d*^ec±(0exp±z(^ae5X^'' + >^%) , (3.4) 

that are complex conjugated to each other for real X and Y, 

c-(0=^, C-{Y\X) = C^{Y\X). (3.5) 

Note also that, for even M, c^{^) contain b^{^) and /^(O as even and odd parts, respec- 
tively, 

c±(0 = &^(0 + /^(0- (3.6) 



12 



As explained in more detail in [4j, the manifest decomposition into positive- and 
negative-frequency parts gives rise to the quantum fields with the creation and annihila- 
tion operators c^(0) that satisfy the commutation relations 

[cH^i),c^m=0, [c-(6),c+(6)] = 5(^1-6)- (3.7) 

In this paper we will be mostly interested in the classical picture, however. 
For the further analysis it is convenient to introduce complex coordinates 

Z^^ = X^^ + ^X^^ = 3fJZ^^ + (3.8) 

The real part ^Z^^ of Z^^ is identified with the coordinates of the generalized space- 
time X^^ that contain in particular Minkowski coordinates. The imaginary part X"^-^ = 
c^2^AB jg required to be positive definite and was treated in [1] as a regulator that makes 
the Gaussian integrals well-defined {i.e., physical quantities are obtained in the limit 
X"^^ — > 0; note, that the complex coordinates Z^^ introduced in [1] are related to Z^^ 
as 

^AB ^ y The space of coordinates Z^^ forms the upper Siegel half-space S^m 

pn [22]. Evidently, -Z G provided that Z^'^ e Sjm and vice versa. 
The variables can also be complexified 

extending the Siegel space to Fock-Siegel space Sjm x C*^. 

The continuations of the functions C^f l3.4p to the Fock-Siegel space i^M x C*^ are 

C+iy\Z) = ^ /t;^^ec+(Oexp(ae^eB^^^) expz(e^3^^), (3.9) 

TT 2 J 

C-iy\Z) = ^ /'dA^ec-(Oexp(-^naeB^^'') exp-^(e^3^^), (3.10) 

TT 2 J 

where c'^iC) are some "Fourier coefficients". Depending on a problem in question, they 
can be chosen to belong to different functional classes. 

The broadest framework is provided by distributions that grow not faster than expo- 
nentially of order two and zero type at infinity, i.e., not faster than exp > 0. In 
this case, c^(0 belong to the space Sy2(^^^) dual to the Gelfand-Shilov 

space ^1/2 (M^O I- 



^ Recall that Sai,...,aM defined in [41] as a space of infinitely difFerentiable functions (t>{xi, ..,xm) 



such that the inequality 



^fei ™/cM , / N 



< CqA\^ . . . A)^' k'^'°" . . . kl^'""" holds for any 



integer nonnegative fc^, and some constants Cq and Ai that depend on (f>. The Fourier dual 

space S'"!' "'""^ consists of infinitely differentiable functions ip{pi, ..,pm) that satisfy the inequality 

Pi^ ■■■ Pm opl^...opif ^^^^ - ^'^^i' ■ ■ ■ ^m' • ■ ■ 1m'°"'' fo'" ^"^y integer nonnegative k„ q, and some 

constants Ck and B, that depend on V- 5i/2(M*^) and 5'i/^(M^^) are shorthand notations for 

M 

Si/2,...,i/2 and 51/2,- ■■1/2^ respectively . 

M 



13 



It can be shown B that C+iy\Z) dSD and C-(3^|Z) (^M) are, respectively, analytic and 
antianalytic in y. The (anti)analyticity of (^C^{y\Z) j C~^{y\Z) in Z & S^m follows 

from the integral representations ^ (I3.10p j ( I3.9p . 

As a subclass, one can require c^(0 to be infinitely different iable functions, that grow 
not faster than exponentially of order two and zero type at infinity. Then, for any Z G Sjm 
the functions c+(0 exp {i HU^bZ^^) and c-(0 exp ( - i h^A^BZ"^^) belong to Si/2{R^^) 
with respect of ^ and hence their Fourier images C^{y\Z) (13. 9p and C~{y\Z) f l3.10p 
belong to ^^/^(IR*^). It is worth to note that, as shown in ^ |13], the class S^/^^M*^) 
plays a distinguished role in the analysis of convergency of power series in the Moyal 
star-product in noncommutative field theory. This is particularly interesting taking into 
account that the interactions of HS fields (for more detail on the role of star-product 
in HS theories see [21 EH and reviews [271 EHl I29j ) is governed by the Moyal star- 
product, which however acts on the noncommutative spinor variables Y"^ rather than on 
the space-time coordinates as in noncommutative field theory. 

Further restrictions may be imposed in the case where the fields have to be normaliz- 
able with respect to one or another norm. This is needed to guarantee that the bilinear 
currents are well-defined. As discussed in more detail in Subsection 14.31 the relevant 
classes of functions c^(0 include Sobolev spaces Ll(M^) and Schwartz space S'(]R*^). 

It is easy to see, that C'^{y\Z) and C^{y\Z) are complex conjugated as a consequence 
of I.e., C+{y\Z) = C-{y\Z) and 

^ ^ +^^ o.,fo.,H V^(y|^) = 0. (3-11) 



dz^^ dy^dy^^ 

^-zh^^]c-(y\Z) = 0. (3.12) 

.dz dy dy J 

The equations (13.111) and (13.121) uplift the massless field equations for (negative)positive 
frequencies to the Fock-Siegel space. The (anti)holomorphy properties of reconstruct 
them in the Fock-Siegel space in terms of theirs boundary values C^(Y\X) at A^Af x IR^^- 
Remarkably, depending on the sign of the second term, the classical unfolded field equa- 
tions (13. lip and (I3.12p distinguish between positive and negative frequencies, the prop- 
erty usually delegated to a quantization prescription. Let us note that this phenomenon 
also takes place in Minkowski setup with (appropriately complexified) unfolded equations 
(I2.22P as well as in the related twistorial world-line particle models [381 EH]. In this case, 
the analogues of the upper and lower Siegel spaces are the forward and backward tubes. 
The fields can be unified into the field 

c{y,y\z,z) = c+{y\z) + c-(y\z) , (3.13) 

that, however, does not possesses definite (anti)holomorphy properties in Z^^. 



*We are grateful to M.A. Soloviev for communicating to us this fact. 
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4 Bilinear currents 



4.1 Current equations 

Let us introduce a conserved current which generahzes that of in a way convenient for 
the further analysis. The key fact is that a differential 2M-form 

w^''\g)= {dWAA(hWBdZ^^-dy'y^'' g{W,y\Z) (4.1) 

is closed in a domain in C^^^^^{Z^^) x M*^(IVs) x C^'^{y^) provided that g{W, y \Z) 
is holomorphic in the variables y and Z and satisfies the following current equations 

( dl^ + ^^(^^yw) ^ 1^) = • (4-2) 

Indeed, from ( lO) and ( lO) it follows that 

= (''ly^glr- - (fi»'B<ii:'"'-d3''')g|j) An3'""'(gay, ^1^)) =0 

because 

dWcA [dWAA (hW Bd Z"^^ - dy^Y)^' = Q 

and 

As a result, on solutions of (14. 2p . the charge 

Q = Q{g) = [ uj^'^g) (4.3) 

is independent of local variations of a 2M-dimensional integration surface T?'^ . In par- 
ticular, for functions that decrease fast enough at space infinity, it is independent of the 
time parameter in M.Mi thus being conserved. 

Since (14. 2 p is a first-order PDE system, the space of its regular solutions forms a 
commutative algebra 7^, i.e., a linear combination of products of any regular solutions 
of (14. 2 p is also a solution. The algebra TZ is formed by functions rj of the form 

r]{W,y\Z)=e{WA, y"" -HZ'^'^Wb) (4.4) 

with arbitrary regular e{W,y). An extension of this property to the space of singular 
solutions S is that S forms an 7^- module, i.e., although it may not be possible to 
multiply singular solutions with themselves, they can by multiplied by regular ones. 
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To make contact with the currents (12.171) note, that Eq. (14.11) gives rise to conserved 
currents for g{W, y\Z) of the form 

g{W, y \Z) = r^iW, y \Z)f{W, y \Z) , (4.5) 

where rjiW, y \Z) (14.41) is a polynomial solution of (14. 2 p identified with a HS symmetry 
parameter and 

f{W,y,\Z) = (27r)-^^/2y ^Mjj exp(-i WcU^) T{U, y \Z) (4.6) 
is a solution built of massless fields via the generalized stress tensor 

T{u, y \z) = c+{y - u\z) c-{u + y\z), (4.7) 

where C'^iyiZ) satisfies (13.111) while C^{y\Z) satisfies (13.121) . (For more accurate defi- 
nition that respects necessary analyticity properties see Subsection 14. 2[ The appropriate 
classes of functions C^{y\Z) and C~{y\Z) will be specified in Subsection 14.31 ) The 
inverse transform is 

T(f/, y \Z) = (27r)-*^/2 j jjc^^ j^yy^ y _ ^4 

The equations (14.21) for f{W, y\Z) translate to the following equations for the stress 
tensor 

i.e., fiW, y\Z) satisfies (14.21) provided that T{U, y\Z) satisfies (14.91) and vice versa. 
One can make sure, that the bilinear tensor (14.71) satisfies (14.91) by virtue of (13. lip . (13.120 . 
Note that, up to a factor of ih, the equations (14.90 at real X = ^Z, Y = ^y 

coincide with the rank— 2 unfolded equations of [10]. In particular, (14.100 implies 

^3 ^3 ^3 Q3 

dU^dU^dX^D + dU^dU^dX^B 'OU^dU^dX^^ ~ dU^dU^dX^B 
which equation was used in [5] to prove that the form (12.171) is closed. 

4.2 Siegel strip and bilinear currents 

To define the integration around singularities in Section [5] via a deformation of the inte- 
gration surface over X"^^ in to the complex space, we now introduce a generalized 



\t{U,Y\X)=Q, (4.11) 
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stress tensor that depends on positive- and negative-frequency solutions C'^{y\Z) and 
C~{y\Z) of the unfolded equations (13.111) and fl3.12p . respectively, and possesses proper 
holomorphy properties in G^^\z) x C^^{y), where a domain C Sjm will be specified 
below. In the rest of this section we set h = 1. 

For a positive definite symmetric matrix we introduce Siegel 7i-strip ©j^'* C Sjm 
as follows : 

Z^^ ^ ( (TC — '^Z)^^ is positive definite 

^ ' 1^ '^Z^^ is positive definite 

which is a generalization of a strip < 532; < 7i in C to the Siegel space S^m- 
Since C^{y\Z) is anti-holomorphic in 6^^'' ^ C'^^ , one can see that 

c:^{y\z) = c-{y\z-tn) (4.12) 

is holomorphic in 6^^'' x C*^ . Note, that is a parameter of and 

lim C^(r|X) = C-{Y\X) V X G A^A/ , Y G M^'^. 
W— >o 

It is easy to see that a generalized stress tensor 

TniU, y\Z) = C^{y - U\Z) C:^{y + U\Z) (4.13) 

solves (14. 9 p and is holomorphic in ©j^'' x C^^^. Up to a linear change of variables, its 
restriction to the real subspace in the limit 7i — >■ gives the generalized stress tensor of 
[5]. To reproduce the charge (I2.19p we proceed as follows. Let the integration surface be 

S?f^ = a*^(X) X MJ'UW) 

{y=yo} 

where {X) is any M— dimensional surface that belongs to the real subspace of 
defined by the equation { '^Z = uTi.}, where < < 1 is a free parameter. Substituting 
the 2M-form w'^^' {g) fHTTl) with the function g fj45l) into fji^ . we obtain 

Qh- j {dWA^WBdx^''Y' 7]{w,yo\x + ivn)f{w,yo\x + ivn), (4.i4) 

whence, using the Fourier transform (14. 6p . we get in the limit 7^ — 

Qo ~ j {dWA^ WEdX^"")'' J d^'U exp {-iWcU^) r^{W, y^ \X)To{U,yo\X) ~ 

(4 

j eA,...^,,dX^^^^A... A dX^-^- (-.A, 3^o I ^ . . . ^ T(^, y,\X) 



u=o 
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For the parameter t] of the form fl4.4p with polynomial e and 3^o = this gives (12.191) . 

Alternatively, one can integrate over d'^'^y at fixed Z. For example let S^^'^ = . 
{Z = Xo, y = Y} for some real Xq. We have in the limit Ti 

Qo ~ j (dWAAdY^)^ j d^U exp {-i Wc U^) e {W, Y^ - Xq^^ Wb) To{U, Y \ Xq) 



R2M 



j eM...A,,dY^^ A---ArfF^^^£(^-^^, +^Xo^^^^ To([/, F | Xo)| . (4.16) 



For C+ and C of the form (13. 9p and (I3.10p . respectively, we have 
To(f/, F|Xo) = 

J d''i d^'X c+(Oc-(A) exp z [x^'^iU^B - XaXb) + Y^'iU - Aa) - f/^(Ac + 
So we obtain from (14.161) 

0„^/.«5."Ac-tt).-(A,.(-A-«.i(|-|)) 

/ ^ 5^^) ^'(« + ") - L.„ - / ^" ^ (-^«- sl) ^"(«) ■ ("^' 

As expected, the result is independent of Xq. Up to a rescahng of arguments it reproduces 
the alternative expression for the charge obtained in [5] . 

The important improvement of the form of the current (14. 3 p compared to (12.170 is that 
it allows us to introduce singularities (fluxes) in the spinning variables PFq,, Wa', 3^" , 3^" , 

^al3^ ^a'f^' ^Y^^^ 

can bring in a proper singularity into the charge (14.150 which, as we show 
in Section [5], is needed to reproduce HS currents of pO] in Minkowski space-time. 

Finally, let us note that the generalized stress tensor can also be constructed from 
arbitrary frequency fields as follows. Let fields Cj{y\Z) of frequencies rij = ±1 satisfy the 
equations 

(5|5B + "6-534^)w=0. (4.18) 

Then the generalized stress tensor 

T{U,y\Z) =Cj(a,{y -U))\Z^ Ck(^(3k{K + y))\Z^ (^{ajf = Vj, {I3kf = -rik 



satisfies (14. 9p . Note however that the charge (14. 3 p built of regular parameters (14.40 and 
fields of equal frequencies, that are holomorphic in the same Fock-Siegel space, vanishes 
because an integration surface can be deformed to the respective infinity in the imaginary 
coordinates Z^^ where the fields Cj{y\Z) vanish. For the fields of opposite frequencies 
such a deformation is not possible, that results in a non-vanishing charge (may be, after 
an appropriate ?i— shift). 
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4.3 Appropriate classes of solutions 

The formula ( ]4.17p for the charge was obtained under assumption that the integrals 
under consideration are well-defined. The following four options are most significant 

(i) c±(0 e L,iR^') , 

(ii) c±(0 e LliR'-^) , 

(iii) c±(Og5(m^0, 

(iv) c±(0 e Sy,{R'-^) and c^(0 G S[/,{R''^) . 

The cases (i)-(iii) are appropriate, respectively, for the cases of e = 1 (electric charge), 
degree q polynomial e and generic polynomial e ^—2,^, f The case (iv) is appro- 
priate to pair a generalized function with a test one. This is relevant to the analysis of 
the composition formulae for "D-functions in Subsection 16.41 As mentioned in Section [3], 
in all these cases C~^{y\Z) (13. 9p is analytic and C~{y\Z) ( IS.lOp is anti-analytic in y and 
Z provided that Z G 

Then the generalized bilinear stress tensor fl4.13p is 

TniU, Y\Z) = C+{Y-U\Z) C-{Y+U\Z - iH) = (4.19) 

Its Fourier transform 04.61) is 

fniW, Y\Z) ~ (4.20) 
/ ^V+(^^^)^'(^^^) e^P^( - ^WaXbZ^''+ h'-{W + x)a{W + x)Bn^'' +XaY^). 

Note, that for all cases (i)-(iv) 

^)c-(^)eS;/,(K"'). (4.21) 

As a result, since Z G &^m\z), fniW, Y \Z) (ICTj) is integrable over R^ {W) \I{Z, Y) G 
&^l^\z) X R*^ . Moreover, the exponential factor with Z G &^l^\z) guarantees that 
fni^: y\Z) is analytic in 3^ and Z. 

As an illustrative example, let us consider functions of the form 

c^(0 = ^^^(0 expi-E^^^U^s) (4.22) 



5 Recall that L«(M^^) is the Sobolev space, and S(R^'^) is the Schwartz space of infinity differentiable 
functions /(^) that decrease at infinity with all their derivatives faster than any multi-degree of jj-j . 
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with some polynomial -P^(0 ^md positive definite symmetric E"^^^ , that allows explicit 
analysis. Evidently, the integrals (14.171) converge absolutely for any polynomial parameter 
e. Evaluating the Gaussian integral, we have 



fniW, y \Z) = (27r)-^^/2 j d^'U exp (-z Wc f/^) (4.23) 
P+ (z^) det -H-iZ + E+) exp ( - \{-iZ + E+)-\b{U^ - 3^^)(f/^-3^^)) 
det-^(H+^Z+E-)exp(-l(H+zZ + E-)-i^5(t/-^+3;^)([/^+3;^)). 

The real parts 3?(?i + iZ + E~) and ^{—iZ + E^) are positive definite provided that Z 
belongs to the Siegel Ti— strip. Since from JF e S^m follows G S)m (see e.g. [22]). 
3ft((7^ + iZ + E-)-i) and 3fJ((-zZ + E+)-i) are also positive definite. Hence />^(>V, y \Z) 
f l4.23p is holomorphic in its arguments for Z G ©^^^ Note that the "additional" analyt- 
icity of fniy^i y\Z) in W takes place for any c^{C) ^ Si^/2^^^)i which is the case for 



5 From AI4 to Minkowski space 
5.1 Idea of construction 

In [Ij, it was conjectured that the charge (12.191) should reproduce the HS charges (I2.20p in 
Minkowski space [20] by an appropriate reduction to C M.^. Since the charge in 
M-i contains four integrations versus three in the Minkowski space, the naive reduction 
with the fourth integration over a cyclic spin variable in gives zero because the cycle is 
contractible. To make the cycle noncontractible, a singularity (flux) should be introduced 
in the spinning space. As we demonstrate now, this can be done using the generalized 
current (14. ip . which result was hard to achieve starting from the original expression (12.170 . 

As explained in Section [21 the embedding of the Ad Minkowski space-time into A^4 is 
conveniently described in the language of two-component complex spinors. In these terms, 
= {W^, ) and = (t/", 1/°" ). It should be stressed that the complex structures 
of the Siegel space and of two-component spinors are different. Correspondingly, since 
both the real and imaginary parts of the complex variables Z^^ are real symmetric 
matrices, in terms of two-component complex spinors we have 

with 

Note that 3?Z"^ = 3?X"^-S5X"^, = 5>X°^ + 3?X"^, etc. Analogously we introduce 
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Let us choose the integration surface in the form 

E« = (T=^(Z"^') X 3^ - 3^0) X a%W) , (5.1) 

where a^(Z"^') is a three-dimensional surface in the complexified Minkowski space ^M^, 
a^[Z'^^, Z"' ^ , 3^) is a one- dimensional cycle in the "spinning" subspace, a^{W) C M^(Vl^) 
is a four-dimensional surface and 3^o is a free parameter. 

An elementary calculation then shows that the pullback of the differential form fl4.ip 
to the integration surface f l5.ll) gives 

■cxj^^^{g)\^, ~ dZ^yA dZ^^'A dZ^^' A dMW,y\Z) A d^W g{W,y\Z)WpWp> , (5.2) 

where 

K{w, y \z) = w^w, z^"" - w^w,, z^''^' - ly^ 3^^ + w^, y^' . (5.3) 

The key observation is that for any 3^0 the function A(VF', 3^ — 3^o |^) solves (14. 2 p and 
is independent of Z""^ . This allows us to use A to introduce a singularity in a way 
independent of the complexified Minkowski coordinates Z'^^ . 

To obtain Minkowski charge we set 

f7A(W^, 3^ \Z) = A-\W, y\Z) g{W, y \Z) . (5.4) 

Since the functions K{W,y\Z) and g{W, 3^1^) solve the equations (14.20 , the same is 
true for g^iWi y 1^) (15.40 away from singularities. From (15. 2p we have 



j w^^\gf,) ~ j dZay A dZ""^' A dZ^^' a"^ Ad^'W g{W,y\Z)Wp Wp> . 



(5.5) 



S8 S8 



The idea is to choose one- dimensional cycle such that — = i dcj), where (p E [0, 27r) to be 
a real coordinate on a^. This can be done as follows 

(t1(Z°^,Z°''3',3^) = (5. 

{Z"^ = p5"^exp(i0), Z"'^' = pS''''^' exp(i0), 3^" = p5"exp(z0), 3^"' = pS"' exp(i0)}, 

where p > is a real parameter and at least some of parameters S"^, S°'' ^\ S*", and 5*^' 
are non-zero. 

For this choice we obtain 

k{W,y\Z)l, = pexp{t<P)PiW), (5.7) 

where 

P{W) = W^W/jS""^ -Wa'Wp'S''''^' -W^S"" + Wa'S''' . (5.8) 



21 



Hence ^ = id(f), and a residue is at p = for P{W) ^ 0. 

The subtlety that the integrand of the right side of (15.51) is not defined at P{W) = 
does not affect the result for any g(W, y \Z) integrable over ]R^(H^) because P(W) cancels 
out in (15.51) and the variety P : P{W) = has measure zero. As a result, we obtain 

Q ~ y zu^^^Qa) = 2m j £ Z^^' Ad^W g{W, \Z'"'')Wp Wp> , (5.9) 

where rf^Z^^' = dZ^^. A rfZ"^' A c/Z^"' . 
For 

giw, y \z) = r^iw, y \z)f{w, y \z) , (5.10) 

where 'q{W, y \Z) is a polynomial solution (14.41) of (14. 2 p while /(VF, y \Z) is the Fourier 
transform (14.61) of the generalized stress tensor T(f/, y\Z) we have 



1/ f^U 



g ~ j d^ z^^' ^d^w r]{w, yo\z^''')f{w, y^\z)WpWp 
^j'^'^^'m^ (^(^^,3^ol^°'^')m3^ol^)) 1^.0' (5-11) 

which is just the anticipated expression for conserved charge in Minkowski space. 

To give precise meaning to this construction it remains to identify ^^(Z"^') C to 
choose appropriate Siegel strip ?i, replacing the generahzed stress tensor T(f/, y\Z) by 
TniU, y \Z) (I4.13P bilinear in the massless fields C^. This is done in the next Subsection. 

5.2 Integration cycle 

Let a positive definite matrix ?i G a/ be chosen in such a way that = m'^' = 0. 
Let = cr|^(^"^ ) be a real three-dimensional surface in a real subspace of "^M^ defined 
by the condition QZ"'^ = z/7i"^ , where < z/ < 1 is a free parameter. Note, that with 
this choice, both QZ°'^ and 7i°^ — QZ"'^ are positive definite for V2^"^ G o"|^. 
Let a^{Z°'^ , 2" ^ , 3^) be chosen in the form (15. 6p such that both 

= (^!3(p5o"^exp(20)) , lyW^"' , 53(p5o°'^' exp( 

and 

^AB_^AB 

are positive definite, which is true for sufficiently small p because this is the case for p = 0. 
For this choice we obtain that x c 6^?^ ^ C^{y). 

Now let T = T-j-ciJJ, y\Z) be the stress tensor ( ]4.13p built of the massless fields 
with the help of a positive definite regulator matrix Ti. Let fniW^ 3^1^) in (I5.10p be 
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the Fourier transform (gl]) of T = TniU, y\Z). If are of the form ([33]), fl3:T0D with 
functions from one of the classes (i)-(iv) of Subsection I4.3[ then the pullback of the 
function g{W, y\Z) (15.101) to the integration surface (15. ip is integrable over M^(iy) 
and the considerations of the previous subsection are true. Therefore (15.111) acquires the 
form 



(7=0 



In the limit 7i — > 0, this gives for rj (14.41) 
Recall that according to (I4.13P 

To(f/, y\x) = c+(3^ - u\x) c~{y + u\x) . 

Up to a rescaling of variables, for 3^o = this gives the conserved charge in Minkowski 
space of [20] for polynomial e{W, y). 

Let us stress that Q (15.121) is 3^o independent, i.e., -t^Q = 0, because the variation 

over yQ is equivalent to a local variation of the integration cycle away from singularities. 
For e = const and T of the form (14.131) this gives the following identity 

which will be used in the further analysis. 

6 Integral evolution formulae 
6.1 P-functions 



D— functions of the massless field equations (II. ip and (II. 2p in Mm were introduced in [4] as 
their singular solutions resulting from the integral representation (13. 4p with = ^iiT~^ 
(in this section we set h = 1), 



V+{Z) = -— I d^'i e^^iiUisZ^''). (6.1^ 



V-{Z) = -V+{-Z) =V+{Z) (6.2) 
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and 



V{Z, Z) = V+{Z) + V'{Z) = V+{Z) - V+{-Z) . 

By construction, the functions V~{Z), V~^{Z) and V{Z,Z) solve the equations of 
motion fll.ip . For Z in the upper Siegel half-space Sjm, (16.11) gives 

V+{Z) = -i7r-^s-\ (6.3) 

where 

= det{-iZ) (6.4) 

defines a multidimensional hyperelliptic surface and s is chosen to be holomorphic for 
Z = X + i'K^ E S)m {i-e., positive definite X"^^) and to be positive real for purely 
imaginary Z, i.e., X = 0. As shown in [1] 



V+(Z] 



— exp 



x-0 " 4 ^|det(X)| 

where Ix is the inertia index of the matrix X^^ , i.e., Ix = n+ — ri- , where n+ and n_ 
are, respectively, the numbers of positive and negative eigenvalues of X^^. 
From and ([63]) it follows that 

V-(Z) = iix'^s'^ , 

where s is complex conjugated to s (16.41) . 

The dependence of the D— functions on y, y is reconstructed via the unfolded equa- 
tions (13.111) and (13.121) . respectively. In particular, 

v+{y\z) = d''i e^^iiUiBZ^'' + y^U). (6.5) 

leading to 

'^^iy\Z) = ^ eM-\ZABy^y'') , ZabZ^^ = 6^ . (6.6) 
Notice that V+{y\Z) (ESD behaves as -2^H5^\y) at Z ^ 0. 



6.2 Evolution formula in A4i 

The obtained results allow us to give precise meaning to the integral formula (ll.7p in 
A^4 by defining the integration measure as corresponding to the electric charge case in 
Minkowski space. This guarantees that the restriction of so defined integral representation 
to the Minkowski space correctly reproduces the dynamics of massless fields. 

Namely, we apply formula (I5.12p with the generalized stress tensor (14.131) at ^ 0, 
where C~{Y\X) is the restriction of C~{y\Z) to the real subspace while C^{y\Z) is 
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replaced by V^i^y — y \Z ~ Z) with y and Z interpreted as parameters. (Recall that if 
2^1, 2^2 e 9)m then Z\ — Z^ G f)M as well, including the case, where either Z\ or Z^ is 
real, belonging to the boundary of ^a/.) Let the resulting conserved charge (15 . 1 21) with 
£ = |, which is a function of y and be denoted as C~{y \Z). We obtain 



c-{y\z) 



dUPdUP' 



j(V+{Yo-U -y\x' -z)c-{Yo + u\x'- 



u=o 



(6.7) 



where Yq is a free parameter, C M^lX) is an Euclidean three-dimensional subspace 
M'^ of Minkowski space, 

a^X) = {T^p.X^I'' = to}, (6.8) 

where a positive-definite matrix T^pi describes the time arrow, the time evolution param- 
eter is t = |7^/j'X"'^ and to is its value associated with the chosen space surface. For 
example, for = dap', X"^' = td"^' + xVf^ , where crf^ are Pauli matrices. 

The key observation is that from (16. 7p it follows by virtue of (15.131) that all derivatives 
of C~{y \Z) with respect to 3^ are related to those of C~{y\X') with respect to y just 
in the same way as C~{y\Z) and C~{y\X') in (16. 7p . i.e., 



d 



d 



d'X'^^' 



dy^ 

j(v+{Y,-u-y\x' -z' ^ 



dU^dUP' 



dyAr 
d 



-c-{y\z) 



-C-{Y^ + U\X'\ 



(6.9) 



u=o 



Let us now show that from (16. 7p it follows that 



Using that d^X^'^' 



C-{Q\Z) 
d^x TT^' 



C-{Q\Z) 



we obtain from (16. 7p 



C-{<d\Z)=i- I rfV (C-(0|X')I'+(0|X'-Z) -C'-(0|X')P+(0|X'-Z) 



(6.10) 



(6.11) 



where we have taken into account that the functions and C satisfy the unfolded 
equations (13.110 and (I3.12p . respectively, together with the facts that 



d 



V+{-U\X - Z] 







u=o 



as a consequence of (16. 6p and 

g2 



d 



^^^^^^jC^{U\X) = ±zC^iU\X) , C^{U\X) = Q-^C^iU\X) . (6.12) 
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To prove (16.101) . it is convenient to use the complex notation for the space coordinates 



X = + ix'^ 



X = x^ — ix"^ 



such that dxi A dx2 A dx^ = ^dx^ Adx Adx . 

The combinations of dual to the coordinates x^, x, x ^ 



(6.13) 



map M.^/Z2 on M^, i.e. fcs, /c and k can take arbitrary values. The leftover ambiguity in 
the integration variables for a given ki is the overall phase factor C,a ~^ expiiy^^Q, 
ip G [0,27r). (Recall that is identified with —C,a-) We set 



expz<; 



' k 



exp — 



k 



The integration measures are related as follows [1] 

dk^ Adk Adk Ad(l)= -8i{^i^i + G^2)'^6 A c^fi A A '^^2 



(6.14) 



(6.15) 



The map f l6.13p . f l6.14p from M^/Z2 associated with the variables C,a to x described 
by the variables ki, is non-degenerate except for the expected singularity at = 0. 
Note that 

' (6.16) 



ei^i + 6^2 = ^kk + kl = ^Jkf + kl + kl . 
For any C+ ([331) and C' (l3J0l) we have using formulae fl6J[3D -( l6:T6ll 

1 



C±(0|t, x) 



(i ^ c (^) exp ±i{tkQ + x k^ + xk + xk) , 



(6.17) 



c^m, x) 



d"^^ (± z /cq) c^(0 6xp ±i{tkQ + x^ks + xk + xk) 



(6.18) 



Using (I6.17P and (I6.18p . we obtain from (16. lip 



C-(0|to,x) 



4v]-6 



(ix'3 A (ix' A dx' / (ip3 A dp A dp A d'ip dk^ A dk Adk A 



— { h^l c (p, ■?/') exp — i(toPo + a;''^P3+a;'p-|-x'p) expi((x''^ — x^)/c3 + (x' — x)/;; + (x'— x)A;) 

64 VPo koj 



dp^ A dp A dp A dip- — c (p, 'i/') exp i(— toPo — a;^p3 — xp — xp) = C (0|to,a;). 



TT" J »Po 
Analogously, from (16. 9p it follows that 



c-A....A,A^\z) 



(6.19) 
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_ d d ^ 

where C {y \Z)a, a = — • • • ^ — C {y \ Z). This is equivalent to 

dyA^ gyA^ 



c-{y\z) _ =c-{y\z) 



Because the dependence on y fully reconstructs the dependence on Z by virtue of the 
unfolded equation fl3.12l) from here it follows that 

c-(y\z) = c-(y\z), (6.20) 

that proves the evolution formula 

C^{y\Z) = \ J d'X''''^^j^(v^{Yo-U-y\X' -Z)C-{Yo + U\X'))\^^^.{Q.21) 
Analogously one can see that 

c\y\z) = \j d'x''''^^j^{v~{Yo + u-y\x' -z)c+(ro-f/|x')) 1^^^.(6.22) 

Note, that from here the usual formulae for spin and spin 1/2 fields follow 
c-io\z)=i^ J rfV (c-io\x')'D+{o\x' -Z) -C-{0\X')V+{0\X' -Z)'^ , 



C,iO\Z) = I d'x' TP^'C-^{Q\X')Vlp,{Q\X'-Z) , 

0-3 

where P+^,(3^|Z) = ^^V+{y\Z), etc. 

6.3 Evolution formula in the Fock space 

Now let us show that the same result can be obtained for any M from (14.161) at e = 
by the integration over the 3^ variables as 

c-(y\z) = -^ j d''y'v+{y-y\z'-z)c-{y\z') (6.23) 

y=Y' 

with any Z'. Indeed, substituting 

c-{y\z) = I d''i c-(0 exp -i{uy^ + UiBZ^"") 
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into fl6.23p and using the expression (16. 5p for we obtain 

^ J d^'r v^Y'-y\z'-z)c-{Y'\z') = 

I d''^c-{Oexp~^{Uy'' + U^BZ^'') = C-(y\Z). 

6.4 Composition properties of the P-functions 

Let us consider 

c+{z\y) = v+{z\y) and c-(z\y) = v-(z-z'\y-y'), (6.24) 

where Z' E Sjm and y E C^^ are free parameters. One can see, that 
C-(z\y) = I d''i c-(0 exp -^(e^eB^^^ + 3^^a), c-(0 = ^ exp i{UiBZ'^'' + y'^U) 
and 

C^{Z\y) = I d^'^c^iO exptiU^nZ^'' + y^U) , c+(0 = ^ • 

Since c~(^) E S1/2 and c~^{^) E S'1/2 for any Z' E S^m, we can pair and C~ (16.241) in 
the evolution formula (I6.22|) to obtain 

v-iy\z) = \ I d'x''''^J^{vH-u-y\x' -z)i?-(f/|x'))|^^^. (6.25) 

Analogously from (I6.2ip it follows that 

v^{y\z) = \j d'x'^^' ^^^(v-{u ~y\x' - z)v^{-u\x'))\^^^. (6.26) 

These formulae express the composition properties of the P-functions. 

6.5 Solutions associated with the P-functions 

Finally let us note that one can use the formula (16.60 to generate new solutions of the HS 
equations as follows. Since V^{y\Z) solves the HS field equations for any 3^, 

-=4= exp ( - -ZAB{y^ + uV^){y^ + uV^] 
ydet(zZ) ^ 4 
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where is any constant vector and m is a constant parameter, also does. Integrating this 
expression at 3^ = with respect to u with some weight p{u), we arrive at the following 
set of solutions of the equation (11. ip 

C{0\Z) = /(v^V'Zab) det-kzZ) , (6.27) 

where / is an arbitrary double differentiable function. One can check directly that C(0|X) 
(Km solves (O). 



7 Higher rank conserved currents 

In [in], we have introduced higher rank fields in A4m- Positive- and negative-frequency 
rank r fields satisfy the unfolded equations of the form 

{a|5i.±'*.gF;4v}c'*(>'|A-) = o. (7,1) 

where k,l = 1 . . .r and riki is a positive definite symmetric form. In [lOj it was also argued 
that rank r fields in Ai m can be interpreted as resulting from the reduction of rank 1 
fields in J^rM to M.m diagonally embedded into J^rM via 

A particular realization of higher rank fields is provided by products of lower rank fields. 
For example, a product of two rank 1 fields 

C{Y^,Y2\X) = C{Yi\X)C{Y2\X) 

gives a rank 2 field. 

Analogously, we introduce a higher rank generalization of the current equation (14.21) 

{d^+ WHAg^]9iW,Y\X) = 0, (7.3) 

where g{W, Y \X) takes values in Mm x W^'^ x W"^'^ with coordinates X^^, Wja, Yj^ , 
A, B = 1, ...M, j = 1, ...r. Again, particular solutions of the higher rank current equation 
are provided by the products of lower rank currents. The higher rank current equation 
allows us to derive multilinear conserved currents. Indeed one can easily see that the 2rM 
differential form 

r:o,rM[g{W,Y\X)^ = (dW^A A [WkBd X^"" - dYk^^J ^ g{W, Y \X) (7.4) 

is closed provided that g{W, Y \X) satisfies (17. 3p . As a result, on solutions of (17.31) the 
charge 

Q'' = ■CU2rM{g) (7.5) 
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is independent of local variations of a 2rM-dimensional surface YP'''^'^ , In particular, it is 
time-independent hence providing a conserved charge. 

The formula (I7.4p gives rise to conserved currents for g{W, Y \X) of the form (14. 5 p 
where the "symmetry parameter" rjiW, Y \X) is of the form 

7^{W, Y \X) = e{W,A, - X""^ Wub) (7.6) 

and f{W, Y \X) is a solution of (17. 3p related to a multilinear "stress tensor" T(f/, Y \X) 
via the Fourier transform 

f{W, Y \X) = (27r)-'"^^/2 I ^rM^ p^.^ JJjc^^ j,^^^ y ^ ^-^ -^^ 



T{U, Y \X) = (27r)-"^^/2 j (f^W exp {i W^c U^^) f{W, Y \X) . (7.8) 
The equations (17.31) are equivalent to the following equations for T 

i.e., f{W, Y \X) satisfies (17. 3p provided that T{U, Y \X) satisfies (17. 9p and vice versa. 

A 2r— multilinear "stress tensor" T(f/, Y \X) can be constructed analogously to the 
bilinear one [5j as follows 

r 

T{U, Y\X) = \[ C+,{Y, - f/,|X) C-,{U, + r,|X), (7.10) 
i=i 

where C^i{Y\X), . . . , C^riY\X) are solutions of positive- or negative-frequency rank 1 
equations. 

The higher rank form (I7.4p can be interpreted as the puUback of the standard form 
(14. ip in M-rM to the diagonal subspace M.m C M.rM (17. 2p . The multilinear stress tensor 
(I7.10p is nothing but the bilinear stress tensor (I2.18P on the rank one solutions in M.rMi 
that result from rank r solutions in M.m provided by the r-linear products of the rank 1 
solutions in M.m- 

The formula (I7.10p gives rise to multilinear currents built of free massless fields. 
Naively, the higher conserved currents constructed from this stress tensor amount to 
algebraic functions of the bilinear currents. This is indeed true if the integration mea- 
sure factorizes into a product of lower-rank measures, that concerns both the integration 
surface and the parameter function ri{W,Y\X) (17. 6p . but may not be true e.g. for non- 
polynomial rj{W^ that contains a singular dependence analogous to (15. 4p in the case 
of Minkowski bilinear current. Note that the resulting nontrivial currents should be nonlo- 
cal from the perspective of usual Minkowski space-time because the charge (I7.5P contains 
2rM integrations instead of 2M in the rank 1 case. Specific examples of higher-rank 
currents will be considered elsewhere. 
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Note that, analogously to the consideration of Sections [3] and H], all higher rank 
constructions can be complexified. Namely, one can consider holomorphic continuation 
C{y^\Z^^) of the fields C{Yjf\X^^), that take values in iOA/(2) x C'^^^(3^). 

Finally let us note that the higher rank system fl7.ip is invariant under the group 
0(r) X R that acts as follows 

yf^yf = exp(0) Tiyf , ^ Z^^ = exp(20) , 

where G M and Tj-' G 0(r) leaves invariant the metric tensor r]jk- This symmetry can be 
used for the derivation of identities between solutions of the equations (17.11) using the fact 
that if two functions Ci{y\Z) and C2{y\Z) satisfy (17.11) and coincide at some Z = Zq, 
Ci{y\ZQ) = C2{y\ZQ), then they coincide for any Z. 

For example, the following identities hold for the D^-functions 

v^{y,\z) . . . v^{yr\z) = exp (rM0) v^{y,\z) . . . v^{yr\z) . (7.11) 

In particular, for the rank 2 case we obtain 

2^'v^{2U\2Z)v^{2y\2Z) = v^{u + y\z)v^{u - y\z) . (7.12) 

Identities (17. lip are continuous analogues of the well-known generalized Riemann identi- 
ties of theta functions. 



8 Riemann theta functions as solutions of higher-spin 
equations 

8.1 Theta functions in the Fock-Siegel space 

A somewhat surprising property of the massless field equations formulated in the Fock- 
Siegel space Sjm x is that Riemann theta functions form their natural solutions. 
Indeed, a general positive- frequency solution (13. 9p periodic under — > y^+n^, G Z^^ 
has the form 

C+{y\Z)= 4expt{hZ^^{27mA)i27mB) + 27mcy^). (8.1) 

With = 1 and h = ^it^^ this formula gives the standard expression for the Riemann 
theta function [22] 

e{y,Z)= J2 exptn{Z^''nAnB + 2nAy^). (8.2) 

Here the complexified space-time coordinates Z^^ identify with the complex period ma- 
trix that defines quasi-periods of 9{y, Z), 

e{y + mZ, Z) = exp {-iTrZ'^^mAmB - 2mm^3^^) ^(3^, Z) , G Z*^ . 
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Also let us note that theta function is Z periodic in the sense 

where Z^ is any real symmetric matrix with integer elements and even diagonal elements. 

The fundamental reason why theta functions solve the HS field equations is that 
both HS theory [H [3] and the theory of theta functions [22] are based on the S'p(2M) 
symmetry and its Weyl-Heisenberg extension which, on the HS gauge theory side, is just 
the HS symmetry. For example, as mentioned in Introduction, for the case of M = 2 the 
sp(4) ~ o(3,2) identifies with the conformal symmetry of massless scalar and spinor in 
three space-time dimensions. In the case of M = 4, the appearance of the sp(8) symmetry, 
that acts on the infinite sets of all bosonic and all fermionic massless fields, is a less trivial 
fact observed originally in [1]. The equation (11.31) . which is the simplest 5p(8) invariant 
unfolded equation, was shown in [3] to describe properly massless fields of all spins in four 
dimensions (a closely related argument was also given in [2]). On the other hand, it is 
well-known that theta functions form a ri^2-niodule, where ri^2 is the Igusa subgroup of 
5p(2M|Z) [22j. 

Thus, the fact that theta functions satisfy the same equations as HS fields is not 
too mysterious once Sp{2M) appeared in the HS theory. Moreover, from the HS theory 
perspective, the Fi 2 symmetry in the theory of theta functions is the leftover symmetry 
of the continuous HS symmetry that leaves invariant a particular solution of the HS field 
equations provided by the theta function up to a phase. This class of solutions may indeed 
play a distinguished role in the HS theory because conserved currents constructed from 
such solutions, turn out to be invariant under Fi^2- 

The roles of the variables Z^^ and in the HS theory and the theory of theta 
functions is to some extent opposite. In the HS theory, Z"^^ are space-time variables 
while the twistor variables play an auxiliary role at least in the conventional field- 
theoretical picture. The indices A = 1,2 ... M for M = 2^ are interpreted as spinorial on 
the HS theory side. In the theory of theta functions, M identifies with genus g, the period 
matrix Z^^ (usually denoted fl"^^ [22]) plays a role of a parameter, while the dependence 
on y^ (usually denoted z^) is of most interest. Note however, that in the nonlinear HS 
theory it was realized since nineties (see [27] and references therein) that the fundamental 
HS dynamics is encoded in terms of the twistor variables y^, while the role of Z^^ is to 
visualize the HS dynamics in terms of local events [1] . 

Many of the well-known properties of theta functions acquire a nice interpretation in 
terms of the HS symmetry (12. 9p of the fundamental unfolded equation (13. lip which is 
a multidimensional analogue of the Schrodinger equation. For example, theta functions 
with characteristics 6^ G R*^, a a G R^^ 

e['^] {y, Z) = exp [tirZ^^aAaB + 2maA3^^ + 2i7raA&^) ^(3^ + Za + b, Z) 
= exp{inZ^^{nA + aA)inB + aB) + 2i7i{nA + aA)iy^ + b^)) 

that also solve (13.110 . result from the action of the HS symmetry (12. 9p with Ja = 2i7raA, 
= and /i = j- on the theta function. 



32 



Consider the set of theta functions with equal characteristics = = or ^, 
Vyl = 1,2...M, 

e[:]{y,Z)= J2 expm{Z^''{nA + aA){nB + aB) + 2{nA + aA){y^ + a^)), (8.3) 
which consists of 2*^ independent functions. It is easy to see that 

e[:]{-y, z) = (-i)S"-2-e[^](3^, z), 

i. e., 2*^^^ functions with an odd number of = | are odd in and 2*'^"^ functions with 
an even number of a^i = | are even in y^. In accordance with the normal relationship 
between spin and statistics, the odd functions describe half-integer spin massless fields 
while the even functions describe integer spin massless fields (in the former case the 
solution has to be multiplied by a Grassmann odd element). In particular, the theta 
function (18.21) is a member of this set with = 0, i.e., it is bosonic. 

The class of solutions C^{y\Z) (18.11) of the unfolded HS field equations in A^m is 
special in the sense that C(F|0) is not a regular function of Y as is usually assumed in 
the unfolded HS analysis, but becomes a distribution at real Z. 

Indeed, one can see that 

^[a](3^, 0)= [-l)T.'L,^a^n^ ^yA ^ ^A _ ^Ay 

This formula means in particular, that 6'[^](3^, Z) develops a singularity at 2 0. 
As such, it is analogous to the D-function of the massless field equations. In fact, it 
is the D-function of HS field equations for solutions with appropriate (anti)periodicity 
conditions in y. In the limit in which the period of y variables tends to infinity, i.e., the 
Fourier series in (18.21) is replaced by the Fourier integral, 9{y, Z) becomes the D-function 
as it is obvious from (16.51) . The counterpart of the evaluation formula (16.231) for periodic 
in ^y function C'^\Z), conjugated to C+(y\Z) (EI]) is 

C-{y\Z)= j d'^Y' e{Y' -y\Z' -Z)C-{Y'\Z') (8.4) 

with any Z'. (The direct proof is analogous to that of Subsection 16.31 ) The evolution 
formula for C+{y, Z) is obtained with the help of e-(y,'z) = e{y,Z). 

Also, let us note that the generalized Riemann identities for theta functions [22] are 
discrete analogues of the identity (17. lip for V^. For example, the generalized Riemann 
identity 

J2 o[o] (23^, ^2Ho] m 2Z) = eiy + u, z)e{u -y,z) 

2aAe(Z/2Z)A-f 

can be derived along the same lines as (I7.12p . 
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8.2 Theta functions and massless fields in Minkowski space 

In the case of M = 2, eq. ( 18.3P gives solutions for massless scalar and spinor in three 
dimensions. In the case M = 4, it describes the superpositions of massless fields of all 
spins in four dimensions. Let us stress that direct identification of theta functions with 
solutions of id massless field equations turns out to be so simple just because, in the 
Sp{8) invariant framework, massless fields of all spins turn out to be involved. This is a 
manifestation of the general feature that linear and nonlinear field equations for massless 
fields of all spins are in a certain sense simpler than those for specific lower spins. 

The reduction of the theta function solutions to a definite spin in the Minkowski 
subspace of A4m is more subtle. To this end, let us first of all make precise the relationship 
between Majorana indices A,B,... and two-component indices a, (3 . . . and a', Let 

= {A^, A^, A^, A'^). We set for two-component vectors A 

= (A^ + iA\ A^ + iA^) , i"' = {A^ - iA^, A^ - iA^) . 
The inverse relations are 

A' = l{A' + A''), A' = l{A' + A''), 
Introducing 

ui = Til - ins , z/2 = - , z/i/ = Ui + in^ , z/2' = n2 + in^ , 
which describe points with integer coordinates on the complex plane, we observe that 

uaA^ = z/„i" + z/^/i"' . 
Also we introduce the complex characteristics 

«! = Oi — ms , ^2 = 02 — ia4 , ar = ai + ia^ , 02' = 02 + ^04 . 
As a result, the theta function with characteristics (18.31) can be rewritten as 

e[i] {y, z)= exp + 2Z"^'/ia/i/3' + ^"'^Va'/i/3' + 2^.3^" + 2/i«'5^"') , 

where 

The solutions of spin s field equations in the complexified Minkowski space, that result 
from 6'[^] (3^, Z), are 

Cai...a2si.Z) = (2z7r)^' ^ /iai • • • fJ'a2s exp (2z7r2:°^'/ia/^/3') , 
Ca'i...a'2siZ) = {2i7rY^ ^ Ha'i ■ ■ ■ l^a'is ^Xp {2mZ''^' flafliJ,) . 



2r ^' 2r ^ 
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9 Conclusion 



The main technical result of this paper is the definition of the proper integration mea- 
sure for the conserved charges of 4d massless fields in terms of an integral in the ten- 
dimensional matrix space Ai^. This allowed us not only to reproduce the known HS 
charges in Minkowski space starting from the ten-dimensional matrix space Ai^ but also 
to give the integral evolution formulae for massless fields via the P-functions in Ai^. The 
precise integration prescription is given in terms of the Siegel upper half-space Sjm [21] 
of complex M x M symmetric matrices = Z^^ with positive definite imaginary 

parts. More generally, we observe that massless fields are most conveniently described in 
terms of the Siegel space with complex matrix coordinates. In this setup, positive- and 
negative-frequency solutions are described, respectively, as holomorphic and antiholomor- 
phic functions in the Siegel upper half-space f)M- The systematic reformulation of the HS 
fields in the Siegel space leads to a number of surprising conclusions. 

One is that the unfolded form of the classical massless field equations studied in this 
paper distinguishes between positive and negative frequencies, i.e., particles and an- 
tiparticles, the property that is usually delegated to the quantization prescription. We 
interpret this intriguing observation as the important indication that the unfolding pro- 
cedure is able to describe quantization. It is worth to note that the unfolded equations 
themselves have a form of a multidimensional Schrodinger equation. 

Another important observation is that Riemann theta functions provide a natural 
class of periodic solutions of the 4(i massless field equations. This fact is a consequence of 
the S'p(8) and HS symmetries of the massless field equations. The setup of unfolded HS 
field equations is convenient for the analysis of properties of theta functions. For example, 
Riemann-type identities [22j can be derived by analysing solutions of the rank r equations 
(17.11) analogously to the analysis of massless field P-functions in Section [61 

Theta functions provide non-zero solutions of massless field equations that are in- 
variant up to a phase factor under the transformations from the Igusa group Fi 2 C 
Sp{2M,Z) [221 . As a result, theta functions may indeed play a distinguished role in the 
HS theory as highly symmetric nontrivial solutions because the observables constructed 
from such solutions, like, e.g., conserved currents, turn out to be invariant under Fi 2. 

An intriguing possibility would be if such a solution can be identified with a nontriv- 
ial vacuum of the HS theory. Although any such a solution breaks down the Lorentz 
invariance to some its discrete subgroup, such a breakdown can be compatible with the 
observations if the scales of the periods of the vacuum solutions are large enough. In that 
case, the breakdown of the Lorentz invariance may have cosmological implications. 

The fact of natural appearance of theta functions in the HS gauge theory is expected 
to shed light on a still mysterious relationship of HS gauge theories with String Theory 
and integrable systems. 
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Appendix. Ring of solutions of massless equations 

As observed in Section HJ the space of solutions of the first-order unfolded current equa- 
tion (14.21) forms a commutative associative algebra. A less trivial fact is that the space 
of solutions of the second-order unfolded equation (14.91) can also be endowed with the 
structure of associative commutative algebra as follows. 

Let us introduce the following commutative and associative product o on the space of 
fields A(Y\X): 

(AoB) {Y\X) = ^(V'|X)exp|-2/i^X^^^| B{Y\X). (A.l) 

The space of solutions of the unfolded system (11.31) is closed under the product o, i.e., 
given two solutions A, B of (11.31) . AoB is its new solution as is easy to see by straight- 
forward substitution of (]A.1|) into (II. 3p . 

The meaning of the product o is quite simple. It corresponds to the usual product of 
the "initial data" A(y|0) and B{Y\Q) of the respective problems as follows from (lA.ip at 
X^^ = 0. So, it is not too surprising that the space of solutions forms such an Abelian 
algebra. It is remarkable, however, that the product o has the simple and constructive 
form flA.ip in Ai^j. The integral version of the formula flA.l|) with /i = 1 reads as 

(AoB) {Y\X) = ^^^^ I dWdVA{Y + W\X)B{Y + V\X)exp-^W^V^XAB , 

where the matrix Xab is inverse to X^^. 

The generalization to rank r fields in Mm [10] is straightforward 



exp 



^dY--' dY^ 



References 

[1] C.Fronsdal, "Massless Particles, Ortosymplectic Symmetry and Another Type of 
Kaluza-Klein Theory"., Preprint UCLA/85/TEP/10, in Essays on Supersymmetry, 
Reidel, 1986 (Mathematical Physics Studies, v. 8). 



36 



[2] I.Bandos, J.Lukierski and D.Sorokin, Phys.Rev. D61 (2000) 045002, 
[hep-th/9904109P . 



M.A.Vasiliev, Phys.Rev. D66 (2002) 066006, [ jhep-th/0106149p . 

M.A. Vasiliev, "Relativity, Causality, Locality, Quantization and Duality in the 



Sp{2M) Invariant Generalized Space-Time", hep-th/0111119 in the Marinov's 
Memorial Volume, M.Olshanetsky and A.Vainshtein Eds, World Scientific, 2002. 
M.A. Vasiliev, Russ. Phys. J. 45 (2002) 670 {Izv. Vuzov, Fizica 45 (2002) N7 23), 
[hep-th/0204167|] . 



I.A.Bandos, Phys. Lett. B558 (2003) 197 | jhep-th/0208110P . 

V.E.Didenko and M.A.Vasiliev, J.Math.Phys. 45 (2004) 197 rhep-th/0301054p . 



M.Plyushchay, D.Sorokin and M.Tsulaia, JHEP 0304 (2003) 013 |hep-th /030T067p ; 

"GL flatness of OSp(l — 2n) and higher spin field theory from dynamics in tensorial 
spaces", hep-th/0310297, 

M.A.Vasiliev, "Higher- Spin Theories and Sp{2M) Invariant Space-Time", 



hep-th /030 1235[ 



O.A.Gelfond and M.A.Vasiliev, Theor. Math. Phys. 145 (2005) 35, |hep-th/0304020 



A.A.Zheltukhin and D.V.Uvarov, Phys. Lett. B565 (2003) 229, [ lhep-th/030415ip . 
I.Bandos, P.Pasti, D.Sorokin and M.Tonin, JHEP 0411:023,2004, 



[hep-th/0407180]. 

I.Bandos, X.Bekaert, J. A. de Azcarraga, D.Sorokin and M.Tsulaia, JHEP 0505:031 



(2005), [ lhip^th/0501113P . 

E.Ivanov and J.Lukierski, Phys. Lett. B624 (2005) 304, qhep-th/0505216P . 
S.Fedoruk and E.Ivanov, Class. Quant. Grav. 23 (2006) 5195, [hep-th/06041in . 
S.Fedoruk, E.Ivanov and J.Lukierski, Phys. Lett. B641 (2006) 226, 
[lhep-th/0606053]. 



P. West, Phys. Lett. B650 (2007) 197, |hep-th/0701026 



E.Ivanov, "Nonlinear Realization in Tensorial Superspaces and Higher Spins", 



hep-th/0703056 



M.A. Vasiliev, Nucl.Phys. B793 (2008) 469. larXiv: 0707. 10851 [hep-th] . 

O. A. Gelfond, E. D. Skvortsov and M. A. Vasiliev, Theor. Math. Phys. 154:294-302 



(2008), qhep-th/0601106p . 

CL. Siegel, " Sympletic Geometry", Academic Press, 1964. 

D. Mumford, "Tata lectures on theta", Boston-Basel-Stuttgart: Birkhauser, vol. I, 
1983, vol. 11, 1984. 



37 



[23] M.A.Vasiliev, Ann. Phys. (N.Y.) 190 (1989) 59. 
[24] M.A. Vasiliev, Mod. Phys. Lett A7 (1992) 3689; 

S.F. Prokushkin and M.A. Vasiliev, Nucl. Phys. B545 (1999) 385, 

r hep-th79806236j . 

[25] M.A. Vasiliev, Phys. Lett. B 243 (1990) 378; Phys. Lett. B 285 (1992) 225. 



M.A.Vasiliev, Phys. Lett. B567 (2003) 139, [hep-th/0304049] . 
[27] M.A.Vasiliev, Contributed article to Golfand's Memorial Volume "Many faces of the 
superworld" , ed. by M.Shifman, World Scientific Publishing Co Pte Ltd, Singapore, 



2000, hep-th/9910096 



[28] A. Sagnotti, E. Sezgin and P. Sundell, "On Higher Spins with a Strong Sp{2, M) Con- 
dition", Proceedings of the First Solvay Workshop on Higher-Spin Gauge Theories 



(Brussels, May 2004), [hep^th/0501156, 

X. Bekaert, S. Cnockaert, C. lazeoUa and M.A.Vasiliev, "Nonlinear Higher Spin The- 
ories in Various Dimensions", Proceedings of the First Solvay Workshop on Higher- 



Spin Gauge Theories (Brussels, May 2004), |hep-th/0503128, 
[30] D. Sullivan, Infinitesimal computations in topology, Publ. Math. IHES 47 (1977) 269 
[31] R. D'Auria and P. Fre, Nucl. Phys. B201 (1982) 101, [Erratum-ibid. B206 (1982) 

496] 

[32] M.A.Vasiliev, Fortsch. Phys. 36 (1988) 33. 
[33] E.S.Fradkin and M.A.Vasiliev, Ann. of Phys. 177 (1987) 63. 
[34] S.E.Konstein and M.A.Vasiliev, Nucl. Phys. B312 (1989) 402. 
[35] E.S.Fradkin and M.A.Vasiliev, Int. J. Mod. Phys. A3 (1988) 2983. 

S.Fedoruk and E.Ivanov, 'New model of Higher- Spin particle", hep-th/0701177 



[37] I. Bars and B. Orcal, Phys. Rev. D75 (2007) 104015, [arXiv/0704 . 0296] . 

[38] T.Shirafuji, Progr. Theor. Phys. 70 (1983) 18. 

[39] I.Bengtsson and M.Cederwall, Nucl. Phys. B302 (1988) 81. 

[40] S.E. Konstein, M.A. Vasiliev and V.N. Zaikin, JHEP 0012:018,2000, 
[ ]hep-th/0010239p . 

[41] I.M. Gelfand and G.E. Shilov, " Generalized Functions", Vol. 2., Acad. Press, New 
York, 1968. 

[42] M. A. Soloviev, Teor. Mat. Fix. 153 (2007) 1351-1363, [arXiv: 0708.0811]. 
[43] M. A. Soloviev, J. Phys A 40 (2007) 14593-14604, [arXiv: 0708.1151]. 



38 



